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HIGH ALTITUDE REMOTE SENSING SURVEYS 
OF LAKE ONTARIO 


(IFYGL PROJECT 1F) 


Introduction 


During the course of the IFYGL, a number of high altitude remote sensing 
surveys were carried out in order to investigate the time and spatial 
variations of large scale surface thermal features on Lake Ontario. 


The principal objective of the aircraft mission was to obtain complete 
coverage of Lake Ontario, with an infrared line scanner, on two or more 
consecutive days. 


In order to provide the synoptic coverage of a large lake such as Lake 
Ontario the infrared line scanner must be flown in an aircraft with a high 
altitude capability. The aircraft used for these experiments was a Falcon 


Fan Jet operated by the Canada Centre for Remote Sensing. 


The sensor package used for the surveys consisted of a Texas Instruments 
RS14 line scanner and a Wild RC10O photographic camera. The RS14 line scanner 


operates in the 8-14 micron band and monitors the thermal emission from the 
| lake surface. The dynamic range of the scanner can be selected and fixed. 


Thus, the infrared imagery has the capability of providing quantitative 


information about the scanned terrain. 


The main function of the photographic imagery was to provide information 
on the lake surface features in the visible wavelengths and examine the 


relationships between these and the infrared imagery. 


A total of three missions were flown during the IFYGL field season. 


‘These occurred on July 7, August 28 and 29, and October 29. All surveys were 


flown at an altitude of 10,000 metres. A series of twenty-one flight lines 


_alternating in a north-south mode was sufficient to provide the desired 
coverage. 


The time taken to cover half of the lake was approximately 3 hours, the 


}remainder could only be surveyed after refueling. Thus, the word synoptic 


.applies only to those features which remain relatively unchanged over a time 
;}period of several hours. Shorter period variations (see Elder and Thomson, 
/1970) were not considered in this experiment. 

sSurvey 1 -— June / 

Desertiptton of the Imagery 


A mosaic prepared for the IR data is shown in Figure 1. The most obvious 
[features are the Niagara plume and the thermal bar. Along the north side of 
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Figure 1. 
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Figure 2, Surface tsotherms prepared from ship data 


the lake the thermal bar extends well into the lake. The southern extent of 
this boundary is in the order of 30 km from the north shore. This particular 
discontinuity was also observed on the photographic imagery where a signifi- 
cant water colour boundary coincided with the thermal boundary. The thermal 
gradients on the south shore of the lake are more compressed and exhibit a 
definite gradation in temperature which shows up aS a series of grey tones of 


the imagery. The mean width of the southern thermal bar is in the order of 
5 kn. 


Another important feature evident on the thermal imagery was a strong 
upwelling (see Bukata and McColl 1973) condition on the north shore which 
extended from the Toronto area all the way through to the Prince Edward County. 


Also of interest is the fact that the western basin (east of an Oshawa- 
Niagara axis) shows a much more disordered surface thermal regime than the 
lake as a whole. While the cause can be attributed to circulation patterns 
in the westerm basin this type of inhomogeneity is not always observed by 
conventional ship data such as shown in Figure 2. 


Surface Isotherms 


Since the RS14 IR line scanner can be operated with fixed dynamic ranges, 
a series of densitometric measurements were made on the data in order to 
exploit the quantitative information contained in the imagery. 


When the gain and D.C. level settings of the scanner are constant over a 
number of flight lines it is possible to compare the scanner film density and 
relate it to the surface water temperature over the same flight lines. 


On the June 7 survey it was possible to compare the film density on the 
last part of the mission which in essence covers the eastern basin of Lake 


Ontario. 


A number of film density measurements were made along the nadir line on 
the section of the flight where the gain and D.C. level settings remained 
constant. Ground truth measurements were used to relate two temperatures to 
two density values and then the remaining density values were converted into 


their equivalent temperatures. 


The ground truth measurements were obtained from a number of sources, the 
principal one being the Atmospheric Environment Service airborne temperature 
survey which occurred on the same day as the overflight. Other water tempera- 
ture data were available: from CCIW ship surveys and the IFYGL meteorological 


buoy network. 


The two base temperatures were taken as 3.0°C and 14.0°C. The lower 
temperature was relatively easy to correlate with the imagery due to the large 
areas of 3 C water present at that time of year. The other base temperature 
was chosen because of a relatively large area of 14 C water which was apparent 


on the airborne temperature survey. 


Figure 3 shows a map of the eastern basin of Lake Ontario with isotherms 


Figure 3. 
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Figure 4. Surface water temperature from Atrborne 
Radtation Thermometer Survey 


obtained from the IR imagery. The main temperature features correspond well 
with the airborne temperature survey (Figure 4) and also the ship survey 
carried out for June 6-8 (Figure 2), confirms the main features. 


The sharp temperature gradients across the north-south portions of the 
thermal bar are well marked. Some evidence of the upwelling situation is 
apparent along the north shore where temperatures close to shore are, in 
general, less than 10 C and increase as one proceeds out from shore. The 
two meteorological buoy stations 10 and 11 give mean temperature reading of 
4.0 C and 9.6 C respectively over the time of the over-flight. 


Careful investigation of the errors involved in making quantitative 
estimates from the RS14 scanner data show that the absolute error, irrespective 
of atmospheric effects, is in the order of # 1.5 C. This rather large error 
places severe limitations on the interpretations made for the quantitative 
data. The important point, however, is that the relative temperature struc- 
ture can be reproduced quite well at least under conditons where a large 
range of temperature is present. 


Meteorological Condittons 


The surface temperature structure as observed during the survey is 
strongly influenced by the meteorological conditions before and during the 
flight. 


By 0300 hrs on June 7, a frontal system had passed through the area 
leaving Lake Ontario under the influence of a strong north to north westerly 
flow. Surface winds over the period of the flight were N-NW at 5-8 mps. This 
N-NW flow had persisted for some 13 hours prior to the flight. Under these 
conditions strong upwelling along the north shore is to be expected. This 
steady wind flow would also cause the extension and compression of the surface 
temperature gradients on the north and south thermal bars respectively. 


Survey 2 - August 28 and 29 
Desertptton of the Imagery 


Figures 5 and 6 show the infrared mosaics prepared from the August 28 and 
29 surveys. These two mosaics show Lake Ontario under summer stratified 
conditions. No large surface temperature gradients, such as appear during the 
period of the thermal bar, are to be expected under these conditions. The 
main features observable under stratified conditions are thermal inputs, from 
rivers or urbanareas, and coastal upwelling. 


The thermal inputs are observable only when the river water or industrial 
plume is at a different temperature from ambient. Coastal upwelling, of 
course, will only occur under certain meteorological conditions. 


The two mosaics shown exhibit all these features. The Niagara plume is 
evident on both days and the influence of the river can be followed along the 
south shore for some distance. Other smaller thermal inputs or plumes can be 
seen, i.e. the warm water can be observed streaming out of Toronto Harbour. 
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Figure 6 


The plume from the Lakeview Hydro plant just west of Toronto can also be 


observed on the August 28 imagery. There is evidence of upwelling along the 
north shore on both days. 


The lack of grey scale contrasts on the imagery is a reflection of the 
small temperature range present during this survey. This fact combined 
with the instrument errors make it wurealistic to attempt the derivation of 
an isotherm map of the surface temperature for this particular survey. It 
is, however, constructive to show a plot of isolines of film density for the 
August 29 data (Figure 7). Without any other temperature data the imagery 
indicates the lack of large surface temperature gradients. The lower density 
values along the north shore are an indication of upwelling conditions. This 
particular condition is confirmed with reference to the airborne radiation 
thermometer and meteorological buoy data taken during the August 29 survey. 
A similar comparison for the August 28 data was not possible due to change in 
the dynamic range of the scanner during the flight. 


The most interesting observations on the photographic imagery are the 
surface streaks indicative of internal wave present on the August 28 and 29 
data. On the August 29 imagery these features (see Figure 8) extend from the 
Toronto area alongthe north shore to Prince Edward County. The wavelength 
range associated with the features is between 300 - 1700 metres which is 
similar to previous observations (Boyce and Thomson, 1972) made from high 
altitude surveys. 


Meteorological Condittons 


A relatively strong frontal system had passed into the area by 2100 hrs 
on August 27. During August 28 the front moved across the lake with surface 
winds along the north shore varying between north to north west at 5-8 mps. 
During August 29 the front continued across the lake and by 2100 hrs was well 


south of the New York coast. 


The prolonged N to NW flow is responsible for the upwelling first 
observed on August 28 and which continued through to August 29. | Ground truth 
‘measurements show mid-lake temperatures in the order of 19 - 20 C on August 28 
with readings of 17.0°C at points on the north shore. 


On August 29 conditions were similar with mid-lake temperatures in the 
range 19 - PEC: However, the upwelling was more pronounced and temperatures 
of 16.5°C were observed at points along the north shore. 


The Niagara plume as observed on the imagery was warmer than the ambient 
lake temperatures. Ground measurements on August 28 yield a value of 22.8 C 
as the highest temperature in the plume. It is interesting to note the 
influence of the warm Niagara plume along the south shore due to the light 
westerly wind flow in this part of the lake during the survey period. 
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Figure 7. Isolines of infrared film data for August 29, 
1972 data 
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Survey 3 - October 29 


Desertption of the Imagery 


The infrared imagery obtained during the mission shows, in general, more 
surface thermal detail than the August survey (see Figure 9). Due to changes 
in gain and D.C. level during the flight it was not possible to make densi- 
tometric measurements for this survey. The dominant feature on the thermal 
imagery is the Niagara plume which appears to influence a large part of the 
western basin. It is interesting to note that while the inner thermal 
boundary coincides with a sharp water colour boundary on the photographic 
imagery the outer thermal feature does not correspond to any feature on the 
visible imagery (see Figure 10). The major shore line turbidity features 
observed from the photographic imagery are shown in figure 10. The direction 
of movement, at the time of the overflight is indicated by the solid arrows. 


The airborne radiation thermometer survey on October 30 shows very cold 
surface temperatures in the western basin with temperatures in the order of 
50 along the Toronto-Hamilton shore line. Since no major weather system 
disturbed the area over the period of these two surveys the surface conditions 
on October 29 would be very similar to those of October 30. The cold 
temperatures along the Toronto area make features like the thermal plume from 
the Lakeview plant stand out very clearly on the infrared imagery. 


Conclusions 


The information contained in the imagery reflects the physical 
characteristics of the lake at three different periods in the annual cycle. 
As far as the infrared imagery is concerned, the optimum periods for applying 
this technique are in early spring, when the lake is starting to stratify, or 
the autumn when overturning or strong upwelling situations are liable to 
occur. With consecutive coverage, the infrared survey provides detailed 
information on the movement of large scale surface features which are con- 
trolled by the dynamics of the lake as a whole. 


Another important aspect of the large scale surveys is the quantitative, 
information provided by the infrared line scanner. It is evident that the 
infrared line scanner provides a more realistic indication of the complexities 
of the surface temperature structure than can be obtained by conventional ship 


surveys. 


The photographic imagery from the large scale surveys has a direct 
application in the observation of sediment transport on a synoptic scale. 
‘This type of information may be of particular value in monitoring the 
characteristics of turbid effluents in problem areas or in areas, such as the 
\upper Great Lakes, where data gathering by conventional techniques is 
(difficult and expensive. 
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BOUNDARY LAYER RESEARCH DIVISION 
ATMOSPHERIC ENVIRONMENT SERVICE 


IFYGL STATUS REPORT - JUNE 1973 


Momentum, Heat and Moisture Transfer in the Atmos heric 


Surface Layer Over Lake Ontario (IFYGL Project 28BL) 


Scientists: G. A. McBean, H. C. Martin, R. J. Polavarapu. 


Setenttfte Program 


The scientific goal of the program was to investigate the transfers of 
momentum, heat and moisture in the atmospheric surface layer over Lake 
Ontario. This includes studies of the transfer mechanisms of the fluxes, 
the variation of the fluxes over 24 to 48 hour periods, the gradients of wind, 
temperature and humidity to 12 m, and the energy balance at the water surface. 
Parameterization of the fluxes in terms of single level observations is also 
being investigated. 


The fluxes were measured directly by two systems: namely, the sonic 
anemometer - resistance thermometer - Lyman-alpha humidiometer; and sonic. 
anemometer - thermistor - refractometer. Both on-line computation and 
Magnetic tape recording for later digital analysis were used. Cup anemometers 
were used to measure the wind profile and thermocouples were used for tempera- 
ture profiles between 1.0 and 12 m. Net radiation and surface water tempera- 
ture were also measured. 


The measurements were made from a bottom-mounted tower installed by CCIW 
off Niagara-on-the-Lake. The electronics and recording equipment were housed 
on a CCIW barge about 500 feet from the tower. 


The project output will be primarily estimates of the fluxes at one 
location at selected times plus information on the transfer processes that 
will aid in deducing the fluxes at other locations and times. In addition 
detailed analyses of the flux-gradient relations, total and turbulent energy 
budgets and the statistics of turbulence will be made. 


Review of Operations 


The spring intensive period was to have been May 1 to 14, 1972. However, 
because of delays in equipment and cable installation the first data wasn't 
collected util May 13. The last day of operations in the spring period was 
May 25. The data was collected with three separate systems. The first was 
wind and temperature (both dry and wet bulb) profiles which were logged on a 
punched paper tape system. The anemometers and thermocouples were mounted at 
ta2,5 4, 8.80d.11 m, The-wind profile system started operating on May 13 
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and was in almost continuous operation until May 25. Examples of the wind and 
temperature profiles output data is given in Appendix l. The temperature 
profile system took much longer to become operational and was not functional 
until May 18 and then ran until May 24. The data for this period is not 
complete in that breaks of several hours occur randomly. The data has been 
analysed in 30 minute periods (32 minute periods for temperature) and the 
total number of data periods is 541 (270.5 hrs) of wind profiles and 235 (125 
hrs) of temperature profiles. 


The second system was the sonic anemometer-thermistor refractometer 
fluxatrons. These measured integrated latent and sensible heat fluxes over 
30 minute averages. The first data was collected May 17 and observations 
continued until May 24. A total of 70 periods (35 hrs) of observation were 
made. These were generally during the daytime. Net radiation, air tempera- 
ture (10m), dewpoint (10m) wind speed (10m) and water surface temperature for 
each period were also observed. An example of the fluxes as measured for 
May 22, is given as Appendix 2. 


The turbulence system used a sonic anemometer-thermometer, platinum 
resistance thermometer and Lyman-alpha humidiometer and the data was recorded 
on FM analog tape. About 12 hours of data was collected during the period May 
18-25 with most being collected on May 19, 20, 22. A tabulation of the data 
collected is given as Appendix 3. The data will be digitized and statistically 
analysed. 


The October intensive period was October 1-14. This fall intensive 
period was also delayed somewhat. The main problems were in the cables 
between the barge and the tower. These cables were laid the first week of 
May and by October had been severly damaged by the continual motion of the 
barge. For this reason the observational program in October was quite cur- 
tailed. The wind profile system was in operation almost continuously from 
October 5-14 but the temperature profile systems were not operated at all. 
The temperature system required 110 VAC to operate the aspirator motors. 
Unfortunately the 110 VAC power cable was broken during the summer and a 
replacement cable was also broken within a day of laying it at the beginning 
of October. The type of wind profile data for October will be as per Appendix 
I and a total of 398 periods (199 hrs) was collected. 


Through a variety of cable improvisations it was possible to operate the 
fluxatron: systems: for a’ total’ of 20 hrs on October 6, 105 11, 12281357, @ine 
main trouble was with the sonic anemometer cable which had some intermittent 
breaks. About 5% hours of turbulence data was collected on FM analog tape 
on October 10, 11, 13. For two of those hours the sonic anemometer horizontal 
wind channels were not working. A summary of the data collected is given in 
Appendix 3. 


Summary of Progress and Future Plans 


The data from the fluxatron systems has all been analysed and papers 
based on it were presented by H. C. Martin at the IFYGL Symposium at the Great 
Lakes Research Conference at Huron, Ohio and at the National Congress of the 
Canadian Meteorological Society, Halifax, May 30 - June 1. Abstract of one 
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presentation is given as Appendix 4. 


The profile data has undergone the initial analysis stage (as seen in 
Appendix I) and will be further analysed as per the scientific program in the 
next few months. The turbulence data has been scrutinized and will be 
digitized and statistically analysed in the next few months. 


Calendar summaries of the data by days are given in Appendix 5. A more 
complete data report of this program's results will be available by October 1. 
The main data of interest to the program scientists are other results of the 
Niagara Bar experiment and nearby tower and shore data. 


T-33 Airplane Investigations of the Boundary Layer (IFYGL Project 15BL) 


Scientists: G. A. McBean, E. G. Morrissey (Small Scale Processes Research 
Division) 


Setenttfte Program 


This is a project to investigate the micro and mesoscale variations with- 
in the boundary layer over Lake Ontario. The goals are to estimate the 
spatial variations of wind and temperature over a scale range of 1 to 100 km 
and to investigate the horizontal and vertical variation of the fluxes of 
momentum, heat and moisture. The T-33 airplane, operated by and in cooperation 
with the National Aeronautical Establishment, Ottawa, was flown in flight 
patterns of long legs at one level both along and cross wind as well as 
shorter legs near one location at a series of levels between 20 and 300 nm. 


Power spectral analyses of Doppler winds, temperature and flux data will 
be carried out to investigate the mesoscale variations. The wind data will 
be corrected for side-slip. The fluxes will be averaged over 1 km legs before 


the spectra are computed. 


The turbulence sensors have been used to calculate the fluxes of momentum 
heat, and moisture. Power spectral analysis of these data will also be done. 


The end products will be a number of estimates of the spatial spectra of 
temperature, wind and the three fluxes and estimates of the fluxes at differ- 
ent locations and elevations. The statistics of the turbulence and their 
variation in the horizontal and vertical will also be investigated. 


Review of Operations 


All the T-33 flights were made in the fall of 1973. A total of 8 flights 
were made but only five will be analysed for this project. The flights were 
according to the flight tracks given in Appendix 6 except for flight 8 which 
was a box pattern over the middle of Lake Ontario and flown December 13. The 
other flights were October 8, 9, 12 and 15. This program is being analysed 
initially in two separate parts. The slower response instrumentation such 


18 


as Doppler winds, mean temperatures etc. are being analysed for mesoscale 
studies by the Small Scale Processes Research Division under E. G. Morrissey 
' The faster response turbulence data is being analysed by G. McBean in the 
Boundary Layer Research Division. 


Summary of Progress and Future Plans 


The turbulence data has all been corrected for aircraft motion and tapes 
of u', v', w', T', q' have been computed for each flight segment. Simultaneous 
to creating the tapes the relevant covariances uw, wI, wq etc. were computed 
by direct cross products. Examples of these outputs and data summary are 
given in Appendix 5. A total of 92 flight segments at altitudes from 25 to 
300m were analysed. Most were at flight level about 150m. The analysis of 
the data will continue with emphasis on turbulence structure and variation in 
vertical and horizontal and examination of details of synoptic events. For 
the latter in particular the surface analyses, IFYGL buoy and shore network 
and radiosonde data for the appropriate days are needed. 
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Appendix 1 


Wind Profiles - (19 May; 1972) In M/S Averaged Over 30 Min 
Previous To The Indicated Times. 


e 


a en ee ee ee be ee A eee 
Height of Sensor 


TIME (EDT) BR i OR A PS) oy EE 
__ rea er ie see nt 
0030 1.14 1.16 1.29 1.30 
0100 84 .90 1.03 | 1.08 
0130 .50 54 65 . 66 .66 
0200 64 72 .78 81 
0230 1.06 1.19 1. 36 1.41 
0300 1203 1.06 1.16 1.34 1837 
0330 74 . 86 90 92 94 
0400 1.95 2.04 Daly 2.20 2.25 
0430 1.51 1.56 1.70 1.70 1.83 
0500 1222 1.26 1.39 1.40 
0530 oF 93 1.02 
0600 2.45 2.54 2.68 2e71 275 
0630 97 1.03 1.13 pats 1.17 
0700 1.46 1.49 1.63 
0730 1.74 1.76 1.99 
0800 he ie 1.88 
0830 1.42 1.49 1.60 
0900 1.32 1.35 1.48 1.48 
0930 1.60 1.64 1.76 1.76 1.83 
1000 2.00 2.07, oP 2.34 2.51 
1030 2.98 3.25 3.36 Bert 4.41 
1100 3.84 4.21 4.40 4.89 5.76 
1130 Ae2] 4.53 4.76 5.27 6.03 
1200 4.62 4.93 5.24 5.81 6.33 
1230 4.99 5.81 
1300 4.47 4.85 5.08 5.64 6.27 
1330 4.65 5.56 
1400 6.10 7.96 
1430 3.91 4.34 4.54 Seo 6.14 
1500 2.86 3.24 3.42 3.90 4.78 
1530 1.78 De10 2.34 278 3.51 
1600 Di pas 3.04 3.45 4.04 
1630 2081 Bal? 3.32 3.74 4.35 
1700 1.61 1.89 1.95 2.20 2.93 
1730 2.48 2.89 sel? 3.60 4.41 
1800 1.95 2.28 2.45 2.85 3.63 
1830 2.04 2.34 2,52 2.91 3.68 
1850 O26 2.59 2.70 2.95 3.78 
2130 1.56 1.66 1.93 2.24 2.46 
2200 1.48 1.59 1.84 2.03 2900 
2230 49 59 61 syn . 80 
2300 . 70 . 80 . 88 . 89 90 
oe ee ee I eS 2 1.76 
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Appendix I (Continued) 


Temperature Profiles - (19 May, 1972). Temperature Differences 
are 1. Caleulated from Lower Hetght Temperature - Upper Height 
Temperature. 2. Averaged Over 32 Min. Following the Specified 
Time. 3. In Degrees Celstus. 


i 


TIME DRY BULB DIFFERENCE WET BULB DIFFERENCE 
(EDT) 1-2M 2-4M 4-8M 8-12M 1-2M 2-4M 4-8M 8-12M 
er eee 

0025 «42 ht Sree! —. 42 75 32 PAT - 34 
0110 24 otk eee | -.43 -67 -24 21 25 
0142 1.30 -.40 ety oP -.44 36 32 35 
0214 47 wat LS: yoy -.44 98 36 29 33 
0246 38 -.44 Ae) -.45 88 33 28 32 
0318 1.03 -.42 =e =e 1.04 - 36 29 34 
0355 78 -.45 ceric 3! aso 81 - 30 28 oe 
0427 41 -.45 eG mr: 61 231 ZI 32 
0459 Zo -.45 —s22 -.43 32 28 25 30 
0520 a4 -.45 woz =-38 43 saz - 30 36 
0552 Sead 3 =o zeuP —.42 = 15 22 29 - 30 
0645 -1.11 -.48 ee -.42 - 30 - 30 - 30 
O717 1.04 =. eer Ay -.40 «23 - 30 ~31 
0812 03 -.06 18 -.04 « 02 14 i3 
1025 te ih =F) =e OZ Oe iS 16 13 12 
1057 25 =29 — 29 1 Uy ot9 Ay), -12 
EIZ9 Soe —.28 iss! aero 19 ae 20 14 
1205 =320 -. 30 =e poh) 16 - 16 12 - 08 
1237 ret cee’ Ge ar P =.30 16 16 12 07 
1400 me -.20 -.47 Seidl Le, 34 45 Pa 5 
1432 —s23 =. 16 se -.40 22 24 27 ~31 
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Appendix 4 


Abstract 


During six days in May and four days in October 1972, net radiation 
(n), sensible (H) and latent (LE) heat flux were measured directly on a 
tower in Lake Ontario. Over one hundred half-hour runs were made. 


In May both H and LE are small, of the order of 2 mW cm-2, and the 
energy budget is overwhelmingly dominated by N with a daytime average 
of over 50 mW cm-2. This incoming energy is used to heat the lake. At 
night all components are close to zero. 


In October H and LE are positive in the unstable conditions with 
average values throughout the twenty-four hour period of 5 mW em~2 and 
11 mW cm-2, respectively. The average rate of cooling of the lake is 
estimated at about 13 mW cm which is higher than other estimates. 


; : — Le 
The evaporation rate in October is 11.1 cm month ', in good 
agreement with other estimates based on indirect measuring techniques. 
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Appendix 6 


NAE-T~33 Flight Summary 


Date Time (E.D.T.) Flight Pattem Comments 
= Se ee a EE Pee od ne ee eae ee ee 
Oct. 4 1250-1410 1 Check-out - not 
complete 
1507-1540 2 
Oct. 8 1012-1110 1 good winds 
1344-1455 2 
Oct. £9 1019-1114 1 good winds 
1300-1400 2 large q_ fluct. 
Oct. 10 1411-1430 Comparison with NCAR Buffalo 
Oct.-12 1250-1350 1 soodadq= luct. 
1557-1703 Z 
Oct, 215 1025-1125 1 fairly light 
1340-1440 Z winds 
Oct. 16 1120-1310 3 Ottawa-Sarnia- 
Ottawa 
Dec. 13 Box Pattern 


In the following data summary (Table 1) the units are msec , 
°C or g/kg. 


TAS - true air speed 


SIGU, SIGV, SIGW, SIGT, SIGH are the r.m.s. values of u,v,w,T, 
and H (humidity) fluctuations. 
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Appendix 6 (Continued) 
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Figure 12. On legs indicated by dashes, no data to be collected. 
Legs B-C, H-D, D-E parallel to wind; A-B, C-D parallel to shore. 
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Figure 13. Leg A-B along wind; C-A perpendicular to wind; A-D 
parallel to shore. Triangular pattern A-B-C to be repeated at 
50", 200",-500",-1000', 8000'. A-D at S00". 
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35 
COMMENTS BY THE U.S. DIRECTOR 


This issue covers IFYGL activities from April 1 to June 30, 1973 (see 
fig. 1). Some reports on events in July and August are included. 


The major emphasis during this quarter was on data management. A major 
accomplishment was achieved by the team from the IFYGL Project Office, the 
Lake Survey Center, and the Center for Experiment Design and Data Analysis 
(CEDDA) , who has developed and tested processing procedures for the Physical 
Data Collection System (PDCS), consisting of buoys, towers, and automatic 
meteorological stations. After an industrious 8 months, production of the 
PDCS Provisional Data Base is beginning, and is expected to proceed according 
to the schedule shown in figure 2. The processing procedures used are 
described in a later section of this Bulletin. 


Considerable progress has been made in processing data from the U.S. ships 
Researcher and Advance II for the Provisional Data Base. Rawinsonde data 
processing by the CEDDA/USAF team is also proceeding successfully. Informa- 
tion received from Brad Bean indicates that the RFF aircraft data processing 
Should be completed during the first quarter of 1974. Jim Wilson's work on 
precipitation, for which he is processing weather radar data in combination 
with other data types, should be finished by June 1974. Water-sample chemical 
analyses have been completed at the EPA Rochester Laboratory and will be 
entered into STORET shortly. 


The above data management efforts are not an end in themselves. The 
obvious reason is to provide users with data. But for us to do so requires 
initiative on the part of those interested in receiving data. You are 
encouraged to file data requests and personally contact the U.S. Data Manager, 
Dave Drury, D2xl, CEDDA, EDS, NOAA, Page Bldg. 2, Washington, D.C. 20235. 
Despite a heavy workload, we will do our best to meet your data needs. 


Attention of all United States and Canadian IFYGL scientists is directed 
to the Second IFYGL Symposium planned with IAGLR at McMaster University, 
August 12 to 14, 1974. Manuscripts must be sent in by May 1974. IAGLR has 
agreed in principle to print a Proceedings in two parts, one of which will 
contain the IFYGL papers. You are encouraged to participate, as indicated 


in the letter reproduced on p. 38. 
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INTERNATIONAL FIELD YEAR FOR THE GREAT LAKES 


STEERING COMMITTEE 
Co-Chairmen 

U.S.- W.J. Drescher 

Cdn.- T.L. Richards 


Members AES File: 5964-31-3 (ACHC) 
Cdn.- J.P. Bruce 
pera Ha ATMOSPHERIC ENVIRONMENT SERVICE 
D.F. Witherspoon 4905 Dufferin Street 
OS ee aubert’  DOWNSVIEW, Ontario M3H 5T4 
D.C. Chandler June 27, 1973 
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IFYGL Principal Investigators 


Reporting of IFYGL Results 


Arrangements have been made with the International 
Association for Great Lakes Research that a full three days of 
scientifie sessions at the next Annual Conference of IAGLR will be 
devoted to papers arising from the IFYGL. There will likely be two 
or three other sessions running concurrently with the special IFYGL 
sessions. 


The conference will be held at McMaster University, 
Hamilton, Ontario, August 12-14, 1974, and will be co-hosted by 
the University and the Canada Centre for Inland Waters. The conference 
will be followed immediately by an event sponsored by S.I.L. (Interna- 
tional Society of Limnology) which is holding its triennial Congress 
in Winnipeg, August 26-30. This evert at CCIW will be a pre-Congress 
special symposium (August 15) on Great Lakes limnology, consisting of 
51x invited papers, and will be open to IAGLR conference participants 
as well as to those participating in SIL pre-Congress Tours. 


You are encouraged to complete your IFYGL projects to a 
report presentation stage for this conference. Manuscripts will be 
required by May, 1974. The IFYGL Joint Management Team will assist the 
Conference Program Committee in selecting papers for presentation. We 
are hoping to arrange with IAGLR for joint publication of a special IFYGL 
volume of papers given at this conference. 


Your cooperation and participation in this important 
evert would be greatly appreciated. 


Wee a y a 
\ (Gr? é.\). 

T.L. Richards BE.” Aubert 
Canadian Co-chairman U.S. Co-chairman 


IFYGL Joint Management Team IFYGL Joint Management Team 


cece. IFYGL Joint Management Team 
IFYGL Steering Committee 


Jointly sponsored by Canadian and U.S. National Committees for International Hydrological Decade 
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U.S. SCIENTIFIC PROGRAM 


Progress on the U.S. IFYGL tasks for the period April 1 through June 30, 


1973, is presented below. Project/Panel status reports follow the task 
reports, 


Tasks 


a 


Phosphorus Release and Uptake by Lake Ontario Sediments 


Principal Investigators: D.E. Armstrong and R.F. Harris - University of 
Wisconsin 


Analytical and experimental goals met following the second sampling 
trip include measurement of total P, total inorganic P, organic P, inter- 
stitial inorganic P, NH4Cl-P, NaOH-P, CDB-P, HC1-P, and the P sorption and 
desorption characteristics of the sediment. The sediment samples were 
obtained on the Researcher between November 6 and 9, 1972, at stations Sas 
83, 92, (3981056625 60, 745, 9343 1102and 302 


Data were collected for comparing squeezing and centrifuge-filtration 
techniques of obtaining interstitial water and their applicability for 
in situ separations of interstitial water for inorganic P analysis. 


Data compilation and evaluation for publication was begun. 


Analysis of the sediments has not presented any major problems. 
However, recent data have indicated a need for more frequent sampling to 
avoid changes occurring from long-term sediment storage. Storage is a 
problem particularly for interstitial inorganic P analysis where immediate 
separation of the interstitial water from the sediment is probably 
necessary. 


Net Radtatton 


Principal Investigator: M.A. Atwater - CEM 


The basic data set for the entire Field Year has been fully coded 
and stored on magnetic tape. Delays in receipt of data have slowed down 
computations of radiation fluxes for the full year for the 30-grid-point 
array, but work is proceeding. Details concerning the computations are 
given in the Interim Report of April 1973, ''The Radiation Budget of Lake 
Ontario Including Cloud Coverage, Preliminary Results," by M.A. Atwater, 
J.T. Ball, and P.S. Brown, Jr. Weekly radiation values in that report 
can be updated for 50 weeks of the Field Year by means of Version A (all 
available documentation). Weekly spatially averaged radiation fluxes are 


40 


presented in figure 3. The weekly Spatially averaged net radiation is 


- shown in figure 4. Figure 5 shows observed and computed solar radiation 


fluxes at Brockport, N.Y. 


The solar radiation model is described in detail in "Numerical 
Computations of the Latitudinal Variation of Solar Radiation for an 
Atmosphere of Varying Opacity," by M.A. Atwater and P.S. Brown, Jr. 
This paper was partly supported by the Atmospheric Sciences Section of 
the National Science Foundation and has been submitted for publication. 


RFF/DC-6 Boundary Layer Fluxes 


Principal Investigator: B.R. Bean - ERL/NOAA 


No report. 


Nitrogen Fixation 


Principal Investigator: R. Burris - University of Wisconsin 


No report. 


Profile Mast and Tower Program 


Principal Investigator: J.A. Businger - University of Washington 


We have analyzed the mean profiles of wind, temperature, and humidity 
for selected periods and are also looking at the eddy fluxes for these 
periods. This has required development of techniques to make corrections 
for any misalignment of the sonic sensors with respect to the mean wind 
direction. For the selected periods, we will first compare the fluxes 
of momentum, heat, and moisture obtained from both the mean profile and 
eddy correlation measurements. Satisfactory agreement between the two 
methods will provide a check on our data reduction and analysis procedures 
We will then begin analyzing the remaining data. 


Status of Lake Ontario Fish Populations 


Principal Investigator: J.F. Carr - Great Lakes Fisheries Laboratory 


No report. 
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Figure 5. Weekly averaged cloud cover and weekly averaged daily 
observed and computed insolation for Brockport, N.Y. 
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10. 


Merk 


Matertal Balance of Lake Ontario 


Principal Investigator: D.J. Casey - EPA 


No report. 


Runoff 
Principal Investigator: L.T. Schutze - U.S. Army Corps of Engineers 


First-cut estimates of monthly runoff from the United States basin 
were derived for February and March 1973. Corresponding runoffs from the 
Canadian basin were not completed. We plan to continue making first-cut 
estimates of monthly runoff from these basins until agreement is reached 
on a method of extrapolating gaged data in each of the two countries. 


No information has been received on discharge measurements at the 
mouth of important tributaries to verify methods of extrapolating gaged 
runoff over ungaged areas. 


Evaporation (Lake-Land) 


Principal Investigator: L.T. Schutze - U.S. Army Corps of Engineers 


No progress made due to manpower shortage. Plans are to continue 
preparing first-cut estimates of monthly evaporation as provisional data 
become available, and to refine estimates when data from other investi- 
gators are received. 


Stimulation Studies and Analyses Associated With the Terrestrial Water 
Balance 


Principal Investigator: B.G. DeCooke - U.S. Army Corps of Engineers 


Activity has not begun. 


Land Precitpitatton Data Analysis 


Principal Investigators: L.T. Schutze and R. Wilshaw - U.S. Army Corps 
of Engineers 


No progress made because of manpower shortage. Study of estimating 
methods will be continued during the next quarter. When a method has 
been selected, precipitation for short periods will be computed and 
compared with results obtained by other investigators. 
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12. Transport Processes Within the Rochester Embayment of Lake Ontario 
Principal Investigators: W.H. Diment, G.F. Bonham-Cartexe and 
J.H. Thomas! - University of Rochester 


Films from four current meters retrieved from the lake last fall 
have been digitized, and the results partly analyzed. Three additional 
meters have been retrieved. 


Last year's results from the Synoptic grid stations (temperature, 
electrical conductivity, beam transmittance, pH, dissolved oxygen, and 
chloride) have been digitized,and analysis has begun. 


Modeling efforts are proceeding in parallel with continued data 
gathering. We have predicted the steady-state wind-driven circulation 
for Lake Ontario as a whole, and the Rochester embayment in particular, 
for winter conditions. We have investigated the sensitivity of the 
same model to changes in the geometry of the lake basin, and to 
spatially variable wind stress. A theory for shallow-lake circulation 
with variable eddy viscosity has been developed. We are working on a 
two-layer, vertically integrated model of Lake Ontario and the Rochester 
embayment, in collaboration with P. Katz and A. Kizlauskas from the 
University of Illinois. 


We have collaborated with R.V. Thomann from Manhattan College, who 
is developing a biological and chemical model for Lake Ontario. We 
have supplied him with our model predictions (complete three-dimensional 
circulation for any wind conditions) on magnetic tape. He will use these 
data as input for his model. 


We are also collaborating with Monroe County Health Department on a 
detailed study of water quality on Ontario Beach and are attempting to 
predict water quality as a function of wind history. Hourly samples 
from two stations off the beach are analyzed biologically and chemically, 
and wind data are recorded every 10 min. 


Two reports on de-icing salts were prepared in connection with 
another project (NOAA Sea Grant GH-106). These may be of interest to 
those studying materials balance in Lake Ontario. 


Without the aid of larger vessels and DECCA navigation, retrieval of 
the current meters is a more difficult operation than anticipated. We 
have had some success with a 27-ft cruiser provided by EPA and used in 
conjunction with our Boston Whaler, but 13 of the 20 meters remained in 
the lake as of July 24, 1973. Funds are short, and, once the size of the 
problem can be more fully assessed, we shall seek additional support. 


} Added as Principal Investigators on this task. 
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Sotl Moisture and Snow Hydrology 


Principal Investigator: W.N. Embree - U.S. Geological Survey 


Work continues on the problem of computing monthly changes in soil 
moisture based on the neutron-log data. An initial determination has 
been made, but the results, and possibly the method of analysis, need 
further study. 


Methods of analysis are being formulated pertaining to the use of 
soil moisture and snow data in predicting runoff. For such prediction, a 
model will be developed that should be improved by the addition of soil- 
moisture and/or snow data for the Black River basin. 


Boundary Layer Structure and Mesoscale Circulation 


Principal Investigator: M.A. Estoque - University of Miami 


Analysis has been temporarily discontinued until rawinsonde, radar, 
aircraft, and surface data from other projects can be obtained. Computer 
programs for integrating a two-dimensional and a three-dimensional model 
are being written. Our effort during the next quarter will be concentrated 
on theoretical modeling. We will test our computer programs at the NCAR 
Computing Facility. 


Mesoscale Simulation Studtes 


Principal Investigator: M.A. Estoque - University of Miami 


No report. 


Lake Level Transfer Across Large Lake 


Principal Investigator: C.B. Felscher - LSC/NOAA 


Several alternate methods were identified for isolating and determin- 
ing the probable effect of wind, barometric pressure, and water and air 
temperatures on the slope of the surface of Lake Ontario. Analysis of 
the feasibility of these methods was begun. 


Nearshore Ice Formation, Growth, and Decay 


Principal Investigator: A. Pavlak - General Electric Company 


Completion date for data reduction has been extended from June 30, 


1973, to July 31, 1973, as a result of illness. Analysis is proceeding 
along the revised schedule. 


18. 


19. 
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Advectton Term - Energy Balance 


Principal Investigator: J. Grumblatt - LSC/NOAA 


No report. 


Occurrence and Transport of Nutrients and Hazardous Polluting 
Substances in the Genesee River Basin 


Principal Investigator: L.J. Hetling - New York State Department 
of Environmental Conservation 


Regular water quality network stations: 


The biweekly stream sampling program is continuing on schedule 
and will continue through October 18. Laboratory analyses of the 
first 10 months of data have been completed and the data are being 
checked for accuracy. 


In April a third set of samples was collected for pesticide 
analysis, and a third set of samples for analysis of mercury, cadmium, 
lead, copper, nickel, manganese, chromium, and fluorides was collected 
in May. Two sets of sediment samples have been taken and will be 
analyzed for phosphorus, iron, magnesium, aluminum, and calcium. 
Sediment, pesticides, and heavy metals will be sampled in August and 
October. 


Stream system nutrient uptake: 


Sampling on three streams (Fish Creek, Holcomb; Spring Brook, 
Lima; and Avon Creek, East Avon) began on May 8, 1973, and has 
continued on a biweekly basis. This routine sampling provides 25 
stream and sediment samples and 3 wastewater samples. Plankton and 
macrophyte samples are taken at selected stations along each stream. 
The fourth stream (Mud Creek, Victor) was dropped because the ratio 
of wastewater flow to streamflow was too low. 


An intensive sampling program to be undertaken in early August 
will provide for daily samples for 5 days on each of the three streams. 
The sampling will include the addition of a number of stations so 
that certain inputs to the stream (particularly tributaries and 
field drains) that are not normally monitored can be adequately 
checked. This will give, for the three streams, a total of 42 stream 
and sediment samples, 3 wastewater samples and several plankton and 
macrophyte samples. These will be analyzed as described earlier. 
After the intensive sampling period, regular biweekly sampling will 
be resumed on August 28, 1973, and continue through October 25, 1973. 
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Boundary Layer Flux Synthesis 


Principal Investigator: J.A. Almazan - CEDDA/NOAA 


The United States meteorological buoy wind data were validated, 
based on a lake-land breeze analysis done by Jason Ching of CEDDA. 
Wind data from the United States and Canadian buoys for July 8 to 21, 
1972, were used to compute diurnal variations of the north-south wind 
component and a similar component normal to the northern and southern 
shores of Lake Ontario. The lake was divided into "northern, middle, 
and southern" regions. The results show that a pronounced lake-land 
circulation was superimposed on a diurnal variation of larger scale. 
Analysis of the temperature field in the lake-land circulation was 
begun. 


Hazardous Material Flow 


Principal Investigator: T. Davies = EPA 


No report. 


Remote Measurement of Chlorophyll with Lidar Fluorescent System 


Principal Investigator: H.H. Kim - NASA 


No report. 


Inflow/Outflow Term - Terrestrial Water Budget 


Principal Investigator: I.M. Korgigian =>U.S. Army Corps of Engineers 


The final report on the Lake Ontario outflow measurements has 
been submitted to the U.S. IFYGL Data Center. The task is completed. 


Use of an Unsteady-State Flow Model To Compute Continuous Flow 


Principal Investigator: I.M. Korgigian - U.S. Army Corps of Engineers 


Work has been delayed due to lack of manpower. Plans are to 
measure discharge Simultaneously at two hydraulic sections on the 
St. Clair River to provide basic data for developing the unsteady- 
State flow model. Its proximity to Detroit makes this river more 
advantageous to measure than the Niagara or St. Lawrence. The results 
will be used to construct the model, which, with account taken of past 


measurements on-the St. Lawrence, will then be adapted for computing 
Lake Ontario's outflow. 
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25. Radiant Power, Temperature, and Water Vapor Profiles Over Lake Ontarto 


Principal Investigator: P.M. Kuhn - ERL/NOAA 


Task completed. 


26. Algal Nutrient Availability and Limitation in Lake Ontario 


Principal Investigators: G.F. beeen N. Sridharan, and W. Cowen - 
University of Wisconsin 


Sampling trips to Lake Ontario tributaries were made on April 6, 
7, and 30, May 1, 27, and 28, and June 15, 16, and 17, 1973. Open-water 
lake samples were collected by CCIW from April 26 to 29, and on June 15 
and were transported to Madison. 


Except for long-term nitrogen and phosphorus availability studies, 
analytical work for the quarter was completed. In nutrient enrichment 
experiments with Lake Ontario water, absorbance, fluorescence, and 
carbon fixation rates were used to measure growth stimulation of the 
natural phytoplankton. Selenastrum growth assays were used to estimate 
percent of phosphorus in particulate and soluble forms in New York 
tributaries. Rainfall and bulk precipitation samples received from 
D.J. Casey, EPA, Rochester, were also subjected to Selenastrum growth 
assays. 


Some of the problems in assaying particulate matter in dilute 
runoff samples have been overcome by use of new methods. New procedures, 
such as high-speed centrifugation and/or ultrafiltration may also be 
required to evaluate phosphorus characteristics of the rainfall samples. 


Future plans include data compilation and evaluation. Very limited 
laboratory work will include rainwater analyses and growth assays. 


27. Wave Studies 


Principal Investigator: P.C. Liu - LSC/NOAA 


No report. 


28. Cloud Climatology 


Principal Investigator: W.A. Lyons - University of Wisconsin 


No report. 


$$ 


; Now affiliated with Texas A & M University. 
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Zooplankton Production in Lake Ontario as Influenced by Envtronmental 
Perturbattons 


Principal Investigator: D.C. McNaught - State University of New York 
at Albany 


No report. 


Change in Lake Storage Term ~ Terrestrial Water Budget 


Principal Investigator: R. Wilshaw - U.S. Army Corps of Engineers 


No progress during this quarter. Plans call for computer program- 
ming and data processing to be completed by the end of 1973. Results 
should be available shortly thereafter. 


Sotl Moisture 
Principal Investigator: L.T. Schutze - U.S. Army Corps of Engineers 
The investigation has not started. 
Lack of manpower and of data from other investigators on precipi- 
tation, evaporation, and runoff continues to delay the start of this 
task. When work is begun, water storage on and below the ground 


surface will be determined as a Single storage factor from the land 
water balance equation. 


Testing of COE (Corps of Engineers) Lake Levels Model 


Principal Investigator: E. Megerian - U.S, Army Corps of Engineers 


No report. 


Nearshore Study of Eastern Lake Ontario 


Principal Investigator: R.B. Moore - State University of New York at 


Oswego 


No report. 


34. 
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Internal Waves - Transects Program - Interpretation of Whole-Basin 
Osetllattions 


Principal Investigator: C.H. Mortimer - University of Wisconsin, 
Milwaukee 


During the period April 1 to June 30, we began to apply some of 
the theories for whole-basin internal oscillations, principally those 
put forth by the Principal Investigator, to the interpretation of 
temperature vs. depth information collected during our three cruises 
aboard the Researcher and Advance II. Isotherm depth fluctuations 
are observed on almost all space and time scales from short, high- 
frequency waves (below the resolution of our instruments and sampling 
program) through whole-basin "'adjustments" lasting several days. 
During the August 7 to 11 cruise, the Lake Ontario region was hit with 
a particularly violent wind storm. Cross sections from this period 
show a considerable amount of activity. Prominent in the midlake 
cross sections are a series of waves or "troughs" propagating regularly 
northward from the downwelling on the south shore. A preliminary 
analysis of these features was included in a presentation at the 1973 
IAGLR Conference in Huron, Ohio. 


Work is progressing on the reduction of the last bits of tempera- 
ture/depth data. Tying loose ends together seems to be more time- 
consuming than ploughing through the main body of the data. It was 
recently discovered that during the July 24 to 28 cruise the tempera- 
ture calibration for the undulating fish had shifted around several 
times. This meant determining new calibrations and rerunning the 
isotherm depth finding programs for about ten transects. It has also 
lately occurred to us that additional information can be squeezed out 
of the undulator records. It may be possible in some instances to map 
the temperature structure much nearer the surface than we are now doing. 
It also may be possible, from inspection of time series plots of 
temperature, to filter out high-frequency internal waves and thereby 
determine mean thermocline depths much more accurately. 


The plan for the next quarter is to (finally) finish up the 
undulator data reduction and to produce a report on all data collected. 
This will involve processing one field tape (four transects) that now 
has a missing (unrecorded) track. The idea is to fill in the missing 
bit so as to restore or preserve the proper parity. In addition we 
intend to press on with our interpretation of the results. As a 
by-product of another research project we now have at our disposal 
a numerical technique to predict the characteristics of low-frequency 
internal waves propagating along a long channel of rather arbitrary 
shape. It shows-the Stokesian edgewaves and topographic Rossby waves 
which were neglected by Mortimer in his treatment of a channel with a 


flat bottom. 
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35. Pontoporeia affinis and Other Benthos in Lake Ontario 


Principal Investigator: S.C. Mosley - University of Michigan 


No report. 


36. Pan Evaporation Project 


Principal Investigator: Teal Nordensone - NWS/NOAA 


Class A evaporation pans were reinstalled on schedule at all 
€vaporation stations. Data have been recorded from April through 
June 1973, after which the stations were discontinued. Dew-point and 
radiation data will be added when they arrive. Computation is 
proceeding of shallow-lake evaporation at the three Canadian stations 
for the winter months. Analysis of the United States data is behind 
schedule because dew-point and radiation data have not been received 
from the collocated IFYGL meteorological stations. 


When the necessary dew-point and radiation data become available, 
shallow-lake evaporation computations will be made and, upon receipt 
of change in energy storage and advected energy data, corrections 
applied to obtain Lake Ontario evaporation estimates. 


An aerodynamic equation will be used to substitute for dew-point 
data for the months of April, May, and June 1973 when the meteorological 
stations were no longer in operation. Development of this equation 
depends upon availability of dew-point data for the period April 1972 
to March 1973. 


37. Stmulation Studies and Other Analyses Associated With U.S. Water 
Movements Projects 


Principal Investigators: J.P. Pandolfo and C.A. Jacobs - CEM 


A report entitled "The Numerical Simulation of the Response of 
Lake Ontario to the Passage of a Typical March Cold Front" has been 
completed and mailed to the IFYGL Project Office. With some minor 
modification to the one-dimensional air/lake boundary layer model, 
the specification of any or all of the atmospheric dependent vari- 
ables -- wind (u, v' components), temperature (Ta), and humidity (q) -- 
as a function of height and time was allowed for the simulation 
experiment. Of primary concern was the detailed time and space 
(vertical) structure of the currents, water temperature, sensible and 
latent heat fluxes, and stresses across the air water interface and 


> 


3 ; : Paka 
C. Hoffeditz, NWS/NOAA, will become Principal Investigator on this task 
assOL JULY 1, 1973, 
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within the water column. The results were compared with a previous 
climatological simulation of the lake (see "Numerical Simulations of 
Lake Ontario with a One-dimensional Air/Lake Model" by J.P. Pandolfo 
and C.A. Jacobs). The report contains much useful information, not 
obtainable from field experiments, for other IFYGL investigators 
concerned with the dynamic, thermodynamic, chemical, and biological 
regimes of the lake during the winter months. The model results 
reveal some of the first physically consistent estimates for March 
of time and space variability, orders of magnitude, and ranges of the 
vertical exchange coefficients for eddy viscosity and conductivity- 
diffusivity, sensible heat flux and stress within the water layer, and 
certain components of the interface heat budget. 


Structure of Turbulence 


Principal Investigator: H.A. Panofsky - Pennsylvania State University 


Analysis was continued of coherence and phase relations among 
wind-speed fluctuations measured at Niagara Bar in August 1972. 
Study of several additional runs showed small coherences even when the 
wind was lined up with the masts. The reason is believed to be that, 
in contrast to earlier runs, the air temperature exceeded the water 
temperature (stable stratification). Phase analysis for another run 
confirmed Taylor's hypothesis extraordinarily well. Additional runs 
will be analyzed with emphasis on zero-Richardson number conditions. 
A hypothesis dealing with diagonal coherence (2 m at one tower and 7 m 
at another) will be tested more fully. 


Continued analysis of buoy wind records failed to show any 
coherence among adjacent buoys. At several buoys, however, interesting 
spectral peaks for speed and direction were found for periods close 
to an hour, with wind speed and direction in quadrature suggesting 
longitudinal vortices. So far, wind directions are not known well 
enough to judge exactly the relation between coherence and angle between 
wind direction and the anemometers. Better wind-direction data will 
be available shortly from CCIW. The buoy analysis will be written up 


shortly. 


Results from this task will be discussed at the IAMAP Assembly, 
Melbourne, Australia. 


Airborne Snow Reconnaissance 
Principal Investigator: E.L. Peck - NWS/NOAA 
The following interim reports were published in May 1973: 


Interim Report # 3, "Soil Moisture Measurements Made at 
Calibration Lines on February 28-March 1, March 9, and 


March 28, 1973." 
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Overlooked in the previous quarterly report was a brief summary of 
plans and current progress provided on March 28, 1973, to B.G. DeCooke, 
Cochairman of the Terrestrial Water Balance Panel. 


Arrangements have been to obtain additional snow data from the 
special collection program conducted by the University of New York in 
Support of the radar snow measurement project. A paper on ''Lake Ontario 
Snowfall Observational Network for Calibrating Radar Measurements" is to 
be presented at an international Symposium on snow and ice to be held in 
Monterey, Calif., from December 2 to 6, 1973. The paper is being pre- 
pared in collaboration with J. Wilson of The Center for Environment and 
Man. 


Further processing and analysis of aircraft data will include the 
mission flown on March 28, 1973. Ground survey data will also be 
incorporated into the analysis to provide final snow-water equivalent 
estimates along the aerial reconnaissance network. An interim report 
will contain a summary of the aerial radiation measurements. Analysis 
results will be given in a final task report. 


Optteal Properties of Lake Ontario 


Principal Investigator: K.R. Piech - Calspan Corporation 


As stated in IFYGL Bulletin No. 4, the objectives of this task are 
(a) measuring the optical properties of Lake Ontario, especially with 
reference to spatial and temporal characteristics and photic zone defini- 
tion, (b) providing inputs to the lake heat budget and to chemical and 
biological studies, and (c) comparing and evaluating three techniques for 
optical turbidity measurements: Secchi disk, irradiance meter/transmis- 
someter, and aerial photographic photometry. 


Efforts to date have consisted of planning and completing data 
collection, and preliminary data reduction (see tables 1 and 2). Prelim- 
inary data reduction is essentially complete, and all data have been 
entered into a format that facilitates temporal, areal, and intersensor 
comparisons required to meet the objectives outlined above. Analyses in 
terms of these objectives will begin during the next reporting period. 


All ship data have been reduced to a graphical format. Preliminary 
analyses indicate reasonable consistency between transmissometer, 
irradiance meter and Secchi disk data, although the usual data scatter 
caused by variations in illumination and surface wave conditions is 
present. 
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Reduction of aerial data necessitates establishing the relation- 
ship between densities on the color film layers and lake reflectance. 
An initial step reduces film densities to relative sensor exposure 
through use of a step wedge. 


In the calibration technique, analyses of scene shadow areas 
are used to establish the relationship between sensor radiance and 
ground reflectance. To put it briefly, a linear relationship exists 
between the radiance just inside and outside a shadow. The parameters 
of this relationship are related to the component of radiance due to 
atmospheric scattering or flare and to the ratio of sunlight to sky- 
light irradiance. Thus, the relationship between sensor radiance and 
scene reflectance can be established by measuring sunlight and shadow 
radiance for a set of shadows and determining the parameters of the 
fit to the data set. The shadow calibration technique requires no 
atmospheric modeling (and consequent model parameter measurement) ; 
rather, the calibration proceeds entirely from the sensor record of 
shadow elements within the scene to determine the atmospheric component 
of radiance. As a result, a calibrated sensor record is obtained in 
which sensor response is directly proportional to scene reflectance. 


A typical sunlight-shadow radiance plot for an end point of one 
of the lake tracks shows that the calibration parameters obtained from 
these end-point analyses appear to be very consistent. The atmospheric 
component of radiance obtained from the calibration parameters is 
usually on the order of 30 to 40 percent of the total radiance from the 
lake. Hence it is very important for this component of sensor radiance 
to be determined as accurately as possible if lake reflectances obtained 
from the aerial photography are to be used for optical characterization. 
Imagery across the lake has been reduced to relative radiance plots. 
Data from these plots, reduced to lake reflectance values, will be 
incorporated into the data matrix format described below. 


To facilitate data comparison and analysis, all surface and aerial 
data are being incorporated into a matrix format that will allow rapid 
recall and display of temporal, areal, and intersensor comparisons. 
Sufficient capacity has been maintained in this formating to allow 
incorporation of other data of interest in optical analyses. 


A more complete quarterly report containing plots of the optical 
data was submitted to the IFYGL Project Office. Scientists having 
need for samples of these data outputs should contact the Principal 
Investigator directly. 
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Storage Term - Energy Balance Program 


Principal Investigator: A.P. Pinsak - LSC/NOAA 


No report. 


Senstble and Latent Heat Flux 


Principal Investigator: A.P. Pinsak - LSC/NOAA 


No report. 


Thermal Characteristics of Lake Ontario and Adveetton Within the Lake 


Principal Investigator: A.P. Pinsak - LSC/NOAA 


No report. 


Oswego Harbor Studies 


Principal Investigator: G.L. Bell - LSC/NOAA 


No report. e 


Mapping of Standing Water and Terrain Conditions With Remote Sensor 
Data 


Principal Investigator: F.C. Polcyn - University of Michigan 


No report. 


Remote Sensing Program for the Determination of Cladophora Distribution 


Principal Investigators: F.C. Polcyn and C.T. Wezernak - University 


of Michigan 


No report. 


Remote Sensing Study of Suspended Inputs Into Lake Ontarto 


Principal Investigators: F.C. Polcyn and C.T. Wezernak - University 


of Michigan 


No report. 
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Island-Land Prectpitation Data Analysis 


Principal Investigator: F.H. Quinn - LSC/NOAA 


Precipitation data have been collected continuously at the six 
Lake Ontario stations, except for some intermittent data caused by 
malfunction of four gages. Data tapes through May 1973 are being 
reduced. Data collection and reduction are on schedule. 


Lake Circulation, Including Internal Waves and Storm Surges 


Principal Investigator: D.B. Rao - University of Wisconsin, Milwaukee 


No report. 


Atmospherte Water Balance 


| Principal Investigator: E.M. Rasmusson - CEDDA/NOAA 


Programming of the budget analyses scheme continued during the 
quarter. The software for fitting the basic measured quantities (u, v, 
T, q) as functions of latitude, longitude, pressure, and time is 
essentially complete. 


Preliminary rawinsonde data are expected to become available 
during the next quarter. Using these data we will do the following: 


(a) Compare observed values at the exact position of 
the observation, i.e., taking into account the drift 
of the balloon, with those obtained from the analyses 
procedure. 


(b) Using values of measured parameters obtained from the 
analyzed fields, compare mean values at the actual 
position of the balloon with those directly over the 
observing stations. 


(c) Compile a statistical summary of the actual position 
of the balloon relative to each station as a function 
of pressure. 


These analyses are expected to provide data to improve the 
analysis scheme and to evaluate the importance of accounting for the 
drift of the balloon. 


Software development during the next quarter will be concentrated 
on completion of the programming of the computations of the derived 
budget quantities from the coefficients of the orthogonal functions 


describing the basic fields. 
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Evaporatton Synthesis 


Principal Investigator: E.M. Rasmusson - CEDDA/NOAA 


Work on this project is now awaiting data. These include (a) 
meteorological data from United States and Canadian buoys, (b) evapora- 
tion estimates from the Terrestrial Water Balance, and Heat Balance 
Panels, and (c) evaporation estimates from the Atmospheric Water 
Balance analyses. 


Ground-Water Flux and Storage 


Principal Investigator: E.C. Rhodehamel - U.S. Geological Survey 


Collection of water level data has continued in order to provide 
information for Studying the hydraulic properties and hydraulic 
nature of lake-front materials. 


The computation procedures for determining the monthly changes 
in groundwater storage (AS) have been refined and a trial run has 
been made. Descriptions of the methods of analysis used in determin- 
ing the groundwater flux and S values are being prepared. Plans call 
for continued processing and computation of data for AS values, and 
preparing data for submission to the IFYGL data bank. 


Spring Algal Blooms 


Principal Investigator: A. Robertson - IFYGL Project Office/NOAA 


Analysis awaits the availability of data. 


Ice Studies for Storage Term - Energy Balance 


Principal Investigator: F.H. Quinn - LSC/NOAA 


Data reduction has been completed. Instrument calibration factors 
have been determined and have been applied to the data, which are now 
being verified. 


Lagrangtan Current Observations 


Principal Investigator: J.H. Saylor - LSC/NOAA 


No report. 
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Ctreulatton of Lake Ontario 


Principal Investigator: J.H. Saylor - LSC/NOAA 


No report. 


Phytoplankton Nutrtent Btoassays in the Great Lakes 
Principal Investigator: C. Schelske - University of Michigan 


No report. 


Runoff Term of Terrestrtal Water Budget 
Principal Investigator: C.K. Schultz - U.S. Geological Survey 


Values of the mean weekly flows for the Field Year have been 
computed and will soon be sent to the IFYGL Data Center. A report 
has been written. This completes work on this task. 


Coastal Chain Program 


Principal Investigator: J.T. Scott - State University of New York 
at Albany 


The basic data report on the United States coastal chain program 
required correction for transport calculations and will be mailed in 
August. Canadian coastal chain data were received from G.T. Csanady 
and processing has begun. This includes computation of baroclinic 
geostrophic velocity and breakdown into along-shore and normal-to-shore 
components of velocity. The result will be a comparison of the five 
coastal chains for cross sections of temperature, along-shore velocity, 
and baroclinic geostrophic velocity for a variety of natural ''events" 
that occurred during the three alerts. We are planning a joint report 
with Csanady on certain aspects of one of the interesting events of 
"quasi-steady" flow that occurred in July 1972. When all five 
coastal chain cross sections are examined simultaneously the flow 
patterns indicate the possible existence of a long baroclinic. ("kKelvan’) 


wave. 


Perhaps the major problem with the coastal chain program is in 
the "filtering" of small-scale perturbations from the "steady" flow. 
We had hoped to have more measurements of continuous data from moored 
current meters, but our own measurements were spotty because we could 
not free the spare current meters from the regular chain program. We 
are obtaining continuous records for several anchored (tower) or 
moored meters from CEDDA and CCIW. This aspect of our task will be 


postponed until early 1974. 
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During the next 6 months emphasis will be given to the ''coherence" 
study, in which we will examine all five coastal chain cross sections 
for specific wind-stress periods. We are preparing the cross-sectional 
diagrams for the entire data set, and for specific "events," each 
lasting for 5 to 10 days. Csanady will attempt to prepare maps of wind 
stress and some estimate of air-water temperature difference for the 
appropriate days. 


We also intend to continue our analysis of transport components 
as described in the IFYGL Technical Plan. The calculations have been 


made, and preliminary comparison of means of actual velocity versus 
baroclinic geostrophic velocity will begin later in 1973. 


Analysis of Phytoplankton Composition and Abundance 


Principal Investigator: E.F. Stoermer - University of Michigan 


No report. 


Clouds, Ice, and Surface Temperature 


Principal Investigator: A.E. Strong - NESS/NOAA 


No report. 


Analysts and Model of the Impact of Discharges From the Niagara and 
Genesee Rivers on Nearshore Biology and Chemistry 


Principal Investigator: R.A. Sweeney - State University of New York 


at Buffalo 


No= report: 


NCAR/DRI - Buffalo Program 


Principal Investigator: J.W. Telford - Desert Research Institute, 


University of Nevada 


No report. 


Mathematical Modeling of Eutrophteatton of Large Lakes 


Principal Investigator: R.V. Thomann - Manhattan College 


The primary effort during the quarter was continued sensitivity 
analysis on the three-layered vertical model (LAKE 1). Preliminary 
verification runs were also made. The two main areas investigated 
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during these runs were the effect of vertical dispersion and the 
effect of phytoplankton settling velocity on the growth of phyto- 
plankton. The vertical dispersion in LAKE 1 is time-variable to 
simulate the development of the thermocline. The effect of vertical 
dispersion was analyzed based on a range of values. This analysis 
revealed the effect mixing has on the magnitude and time location of 
the pulses in phytoplankton growth. The sensitivity of phytoplankton 
growth to phytoplankton settling velocities was also investigated. 
The effects of two different settling rates were determined. 


Plans for the next quarter are to continue the sensitivity analysis 
and verification on LAKE 1. 
65. Cladophora Nutrtent Btoassay 


Principal Investigators: G.F. Lee” and W. Cowen - University of 
Wisconsin 


No report. 


66. Sediment Oxygen Demand 


Principal Investigator: N.A. Thomas - EPA 


No report. 


67. Matn Lake Macrobenthos 


Principal Investigator: N.A. Thomas - EPA 


No report. 


68. Exploration of Halogenated Hazardous Chemicals in Lake Ontarto 


yes is 
Principal Investigators: G.F. Lee and C.L. Haile - Texas A&M 


University 


The sampling program was continued and expanded, and groundwork 
was laid for intensive gas chromatographic (GC) analysis of all 
samples collected and "cleaned up." Water samples were obtained from 
the Don, Niagara, Genesee, Oswego, and Black Rivers (major rivers 
feeding into Lake Ontario) in an attempt to add information concerning 
sources of halogenated organics to our ultimate overall data for the 


Now affiliated with Texas A & M University. 


° Added as Principal Investigator on this task. 
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lake. The 101 samples were collected at a depth of 0 to 1 m within 

1/4 mi from the river mouths. At the sites, the samples were extracted 
by upflow (syphon) passage through columns of six DC-200 coated 
polyurethane foam plugs. A flow rate of approximately 250 ml/min was 
achieved. The columns were returned to Madison, where the plugs were 
removed, exhaustively extracted with hexane-ether azeotrope, and the 
extracts were then cleaned up. 


Benthic fauna samples collected in January 1973 by S.C. Mosley, 
University of Michigan, on the Linmmos were extracted and the extracts 
cleaned up for GC analysis. 


Progress was seen in determining optimum gas chromatographic 
columns and column conditions for the bulk of the GC analyses to 
follow. Several columns were prepared, using liquid gases, such as 
OV-1, DC-200, OV-17, and mixed liquid phases, such as OV-17/QF-1, 
DC-200/QF-1, and OV-1/QF-1 coated on 80/100 mesh Gas-Chrom Q, 100/120 
mesh Gas-Chrom Q, or 100/120 mesh chromasorb G packed in glass columns 
5 to 7 ft long with an inner diameter of 2mm. The columns were 
evaluated by use of standard pesticide solutions, pesticide cocktail 
solutions, and some Lake Ontario fish extracts. Columns showing high 
potential for separation with reasonable retention times include those 
prepared with DC-200, OV-17, OV-17/QF-1, and OV-1/QF-1. Column tempera- 
tures of 160-200 °C were used with carrier gas flows of 7-35 ml/min. 


Major among problems Slowing our progress of this quarter has 
been a greater than anticipated difficulty in determining optimum GC 
columns and column conditions for the bulk of the analyses. Although 
this has caused some delay in the analysis, problems have been solved 
to a workable degree. 


Final sampling should be completed next quarter with fish samples 
from major rivers feeding the lake added to provide information concern- 
ing halogenated compound sources for Lake Ontario. These fish samples 
will be extracted and the extracts cleaned up in the same way as the 
lake fish already collected. The bulk of the routine and exploratory 
GC analysis should be accomplished. Preliminary work on several 
confirmation techniques such as GC/MS, perchlorination, and TLC will 
be initiated. 


Bastn Precipitation - Land and Lake 


Principal Investigator: J.W. Wilson - CEM 


Derivation of precipitation data from the United States radar 
program is essentially on schedule. Reduction of the Canadian radar 
data, currently lagging, is expected to be on schedule by the end 
of the summer. The United States Data Sets 1 and 2 are now complete, 
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in terms of editing, error checking, and compacting of the raw radar 


data, and derivation of hourly precipitation totals. There are 


1,150 hours that were not successfully collected on magnetic tape when 
precipitation was occurring over the lake. Of these, 925 hours can be 
extracted from PPI photographs or from overlapping coverage by the 
Woodbridge radar. Based on experimental results in extracting the 
missing data from film, a proposal is nearing completion to enlarge the 
scope of reducing these photographs. 


Comparisons of radar and gage data during snow storms indicate 
the importance of studies to compare Oswego snow network data and 
radar data before establishing procedures to measure the water from 
falling snow. Such studies are being included in the proposal to 
widen the scope of this task. 


Preliminary procedures were established by D.M. Pollock and 
J.M. Wilson for integrating the radar and gage data from all data 
sources to obtain daily maps of the precipitation over the lake and 
watershed. Results of radar and gage comparisons for Hurricane Agnes 
were examined in considerable detail. A report is being prepared. 


Evaluation of ERTS Data for Certain Hydrological Uses 
Principal Investigators: D.R. Wiesnet and D.F. McGinnis - NESS/NOAA 


Melting ice has been detected in Lake Erie by comparing visible 
and near-IR differential reflectance and is being confirmed by 
meteorological ground truth data and concurrent NOAA-2 thermal IR data. 


Using near-synchronous ERTS-1 data, we have determined the 
approximate ground resolution of NOAA-2's Very High Resolution Radio- 
meter (VHRR) by comparing the smallest identifiable ice crack in the 
VHRR image. Cracks as small as 600 m are being detected. They are 
a high-contrast linear type of target. 


Photographic interpretation and analysis of a single Lake Erie 
image (e.g., for February 18, 1973) can, under certain conditions, 
provide meaningful data on ice dynamics. Preliminary estimates of 
ice movement range from 800 m/hr to 300 m/hratd.43. toc0.16tkn) eto the 
west. A.E. Strong (personal communication), using the side lap on 
successive days of ERTS-1 passes, has been able to plot ice movements 
that range from 1500 m/hr to 500 m/hr (0.81 to 0.27 kn). Considering 
that time of onset of motion is not precisely known, these figures 
are in rather good agreement. It is significant that the wind 
inducing this movement was relatively light (10 kn or 18,500 m/hr). 


ERTS-1 imagery is outstanding for ice-pack monitoring of the Great 
Lakes. Frequency of sampling interval (18 days) is, however, a severe 
limitation. 
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Distrtbution, Abundance, and Composition of Invertebrate Fish Forage 
Mechanisms in Lake Ontario 


Principal Investigator: J.F. Carr - Great Lakes Fisheries Laboratory 


No-report« 


Coastal Ctreulatton in the Great Lakes 


Principal Investigator: G.T. Csanady - Woods Hole Oceanographic 


Institution 


Research objectives for this quarter were (a) further consolidation 
of the theoretical framework of lake circulation dynamics and (b) 
continuation of the analysis of the summer 1972 alert period. 


In interpreting water movement observations, specification of wind 
stress from surface meteorological data remains a key problem. One 
theoretical study attempting to throw light on this problem is W.H.O.1. 
Contribution No. 3011, referred to in IFYGL Bulletin No. 7. This 
paper has been revised and accepted for publication in Boundary Layer 
Meteorology. A second study in the same general area has been prepared, 


Further work was done on the summer data. One puzzling feature 
of these is that considerable barotropic velocities exist perpendicular 
to shore. It appears that the coastal jet has a wavy or meandering 
appearance even along a relatively straight shore. Preliminary calcula- 


for Gulf-Stream meanders. Owing to the lack of alongshore resolution of 
coastal observations, however, it is not likely that this problem may 

be conclusively disposed of. A second generation of coastal chain 
Studies should focus on exploring a piece of Straight coast perhaps 

50 km long, at several locations to describe alongshore variations 

of coastal jets. 


Another aspect of sumer; -as-well Tasttheeralt, datatre the peculiar 
Shape of the thermocline when it meets a Sloping bottom, especially 
on occasions of developing upwelling. A nonlinear theory of this 
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Lake Water Charactertstics 


Principal Investigator: A.P. Pinsak - LSC/NOAA 


No report. 


Snow Observation Network 


Principal Investigator: R.B. Sykes - State University of New York, 
Oswego 


No report. 


Lake Ctreulatton Model 
Principal Investigator: J.R. Bennett - IFYGL Project Office/NOAA 


Programming, checkout, and testing of model options is continuing. 


Lake Ontario Invertebrate Fauna List 
Principal Investigator: A. Robertson - IFYGL Project Office/NOAA 


Work is continuing on correction of invertebrate fauna list. 


Natural Distribution and Vartability of Phystcal Proper teas 


Principal Investigators: E. Aubert, J. Harrison, and R. Pickett - 
IFYGL Project Office/NOAA 


The objective of this task is to examine the temporal and 
spatial variability of physical limnological, meteorological, and 
hydrologic variables and derived variables, e.g., energy, transport, 
etc. Techniques will also be developed to produce data required 
for other tasks. An atlas of physical lake properties will be 
prepared. This task will start with July 1972 data. Techniques for 
application to other months will be developed. 


nl 
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This is a new task initiated in May 1973. 
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Projyect Areas 


Boundary Layer - J.Z. Holland, U.S. Panel Cochatrman 


The Boundary Layer Panel held a meeting in Washington, D.C. on June 21 
and 22, 1973. Before the meeting, principal investigators submitted a ; 
summary of their respective data collection dates, types of data collected, 
and project status. Four Canadian and two United States panel members 
attended. Representatives of the data management offices of both countries 
and of the U.S. IFYGL Project Office were also present. The data management 
representatives presented a summary of the data processing and availability 
Status. Each of the principal investigators presented a report on the 
results of his analysis. 
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OPERATIONS AND DATA ACQUISITION SYSTEMS 


U.S. Field Headquarters 


The field operations under the responsibility of the U.S. IFYGL Field 
Headquarters were completed June 30, 1973. All land sites and buildings 
have been cleared of gear and returned to their original condition. The 
six Physical Data Collection System (PDCS) stations were disassembled. The 
equipment was checked and shipped to the Lake Survey Center (LSC) and the 
National Data Buoy Center. Post-calibration activities included temperature, 
dew-point, pressure, wind speed and direction, and electronic measurement 
unit tests at the Rochester facilities. Checks of the radiation sensors were 
completed by the Atmospheric Environment Service, Downsview, Ontario. Pre- 
liminary transfer of all data on cassette tapes to seven-track tapes was 
completed in Rochester. Equipment used in this transfer was then sent to 
CEDDA, Washington, D.C., for further processing of the PDCS data. After work 
at CEDDA has been completed, the equipment will go to LSC. 


Positions of U.S. IFYGL Ships 


Errata -- The station locations listed in table 2 in IFYGL 
Bulletin No. 7 (pp. 66 through 68) were in error in two 
instances. Latitude for IFYGL station identifier 66 should 
be 43°S54'S4"N. This position was occupied on and after 
June 13, 1972. Before that date, position was 43 54'36''N, 
77°30'00"W. Latitude for IFYGL station identifier 70 should 
have been listed as 43°29'24"'N. 


70 


DATA MANAGEMENT 


Physical Data Collection System (PDCS) 


Erratum -- Figure 4 in IFYGL Bulletin No. 7 (p. 71) shows 
that data were collected from IFYGL station No. 27 (shal- 
low tower off Rochester) for February and March 1972. 
This is in error. No data exist for this period from 
Station 27 as the tower was removed from Lake Ontario in 
mid-November 1972. 


In cooperation with the working team representing CEDDA and the IFYGL 
Project Office, the computer facility personnel at LSC have designed and 
written the programs necessary for processing the PDCS data. These programs 
have been tested and are becoming operational. 


The working team at CEDDA has completed the interpretation and 
application of electronics and sensor calibrations and formalized the PDCS 
data conversion process. This process (see fig. 6) starts with 6-min raw 
data in the following forms: (a) real-time data (2-month disks), (b) RCC 
weekly tapes (seven-track magnetic tapes), and (c) cassette backup data. 
After the raw values have been handled and merged through software routines, 
corrections are applied for the various electronic units and sensor posi- 
tions. The corrected raw data are then stored in the Provisional Data Base, 
with a source flag appended as a first digit (0 for real-time, 1 for RCC, a 
d and 2 for cassette data). Data that cannot be automatically processed 
are left out for possible later manual workup. Next, the data are converted 
into scientific units, and the resulting 6-min scientific values are then 
Stored in the Provisional Data Base on disks, each containing 1 month of 
data. These can be written out as users' files at any time. 


The reason for the order of the above process is that the PDCS elec- 
tronics and the 15 sensor types were calibration checked independently in 
terms of input and output. The checks vary from simple "go"/'no go" 
operational checks to well-controlled multipoint calibration. Table 3 
shows the availability of calibration data. 


A first-order correction for particular electronics was applied to the 
raw data by means of the equation 


R' = aR + 8, (1) 


where R is the raw data, R' the corrected raw data, and a and 8 are the 
linear transformation coefficients. Note that R is always a three-digit 
number, 000 to 999, and R' is always expressed as a number from 000.0 to 
999.9. The coefficients a and 8 are functions of time, a being one and B 
being zero for ideal conditions, and are based on the performance checks 
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Table 3. PDCS caltbration data 


eS Re eee eee 
No. calibration check 
01 Air temperature x x x 
02 Air pressure a . x 
03 Pan evaporation - ss 
04 Precipitation x x 
05 Longwave radiation x Xx 
06 Shortwave radiation x x x 
07 Dew point x xX 
08 Wind direction (buoys) x *® ne 
09 Wind direction 
(except buoys) : is 
10 Wind speed - % he 
gt Water temperature x X X 
12 Water temperature 
(evaporation pan) x x x 
1S Current direction (buoys) = “ apie 
14 Current speed (buoys) x Xx <* 
15 Current speed (towers) x x sls? 


x - Data available and used in the PDCS conversion process, 


* - No electronics internal calibration or performance check data 
available; analog-to-digital converter calibration check used in 
the PDCS conversion process, 


** - Sensor calibration data available, but not applied in the PDCS 
conversion process. 


Blank - No data available. 


V2 


(electronic calibrations) and the internal calibrations that occurred 
automatically every 10 hours. Performance checks on all electronics were 
made both before deployment and after retrieval, and in some cases in the 
middle of the Field Year as well. They consisted basically of laboratory 
multipoint calibrations of the signal-conditioning circuitry and the analog- 
to-digital (A/D) converter and an operational check of the internal cali- 
bration feature, the digital circuitry, the recorder, and the radio. Since 
the electronics design relationships were known, the measured and calculated 
outputs were compared. The internal calibrations were two-point checks of 
the signal-conditioning circuitry and the A/D converter. The inputs were 
values simulating the sensors at approximately 90 percent and 10 percent of 
full scale; the outputs were the raw data values, called Rg (t) and Ryo (t). 


Data from the performance checks are tabulated as MD (measured data) 
versus CD (calculated data). This multipoint data set is fitted to a 
straight line by 


CD = yMD + 6. (2) 


This equation is used in transforming the internal calibration values at the 
time of the performance checks to corrected base internal calibration values, 
called Rgp and Rjg. The actual internal calibration values during the 
operation of the electronics are used with the corrected base values (averaged 
from all performance checks) to form coefficients a and 8 in eq. (1). These 
coefficients for dew-point and air-temperature channels on land stations 

also contained linear corrections for differences in circuit design from the 
other stations. This was done so that in the next stage of the data conver- 
sion all stations could be treated alike.) Coefficients a and 8 were arrived 
at as follows: 


Rog ~ Rio = 


~ Rog (t) - Ryo lt) * 


_ Roo lt) Rig - Riot) Roo 


(4) 
Rog (t) - Ry Q(t) 


Based on a selected integer value with the lowest variance for that period, 
the actual internal calibration values, generally representing 1 month of 
operation, were chosen to be used with the corrected base internal calibra- 
tion values. The errors resulting from this technique will be subject to 
study as part of the PDCS final data reduction. 
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The a, 8 correction of raw data consists essentially of starting 
with a linear transformation based on a multipoint laboratory check 
followed ky monthly adjustment of this transformation using a two-point 
automatic check. The assumption is made that the two points from the 
automatic internal calibration at the time of the performance checks 
fall on the best-fit straight line of the performance check data and 
represent a base value from which drift or offset during actual operation 
can be calculated. If the internal calibration mechanism works properly, 
then the linear transformation can be adjusted based on the internal 
calibration values to a point where the transformation becomes the post- 
retrieval best-fit straight line. Note that some sensors did not have 
Signal-conditioning circuitry (nor an internal calibration mechanism). 

In these cases a multipoint check of the A/D converter is used as the 
basis of a and 8. (All data channels used the A/D converter.) 


In conversion to scientific units, data were checked for values that 
were either too high or too low as a result of noise or serious sensor 
malfunction. The scientific data in the Provisional Data Base are indicated 
as follows: 


-999.000 = Internal calibration data 

-999.000 = Missing data 

-888.000 = Data out of range of sensor equation 

-777.000 = Data fail high/low check 

-666.000 = Transmission error; message not decoded properly 
SFXXX.XXX = Scientific data; s = sign, F = flag, X = digit 


The first step in converting the corrected raw data to uncorrected 
Scientific units is the sensor-dependent linear transformation 


Eo ae Sak! asa (5) 


where Eg is uncorrected scientific data and Sg and S,} are the transfer 
coefficients. The subsequent steps, if any, are different for each sensor 


The data in corrected scientific units, called Ej, are a sensor-dependent 
function of the data in uncorrected scientific units as in 


EL = f(E,) . (6) 


15 


An auxiliary step in conversion to scientific units is correction for 
the geographic location or position of the sensor. This includes correct- 
ing all direction readings so that they are referenced to true North, 
modulo 360, and correcting dew-point readings for the lead resistances. 
Atmospheric pressure is not corrected for elevation; it is reported as 
station pressure. Also, 950 mb is subtracted from all atmospheric pressure 
values. Rounding is to the decimal place at least an order of magnitude 
better than the system error. Table 4 summarizes the rounding, resolution, 
accuracy, and ranges for the 15 sensor types. 


The complete Provisional Data Base includes all available real-time 
received data, input from RCC weekly tapes and cassettes, as well as data 
worked up manually to fill in gaps left by the automatic processing. For 
the Final Edited Data Base, after analysis and user feedback, as much erro- 
neous data as possible will be filtered out and tapes generated in a format 
suitable for archiving, including hourly averages. 


The Provisional Data Base will be complete on or about June 15, 1974; 
the Final Data Base, on or about August 1, 1974. Note that monthly segments 
of these data will be available in the interim. 
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